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Abstract. Energetic particle transport in a finite amplitude magnetosonic and 
Alfvenic turbulence is considered using Monte Carlo particle simulations, which 
involve an integration of particle equations of motion. We show that in a low-/3 
plasma cosmic ray acceleration can be the most important damping process for 
magnetosonic waves. Assuming such conditions we derive the momentum diffusion 
coefficient D p for relativistic particles in the presence of anisotropic finite-amplitude 
turbulent wave fields, for flat and Kolmogorov-type turbulence spectra. We confirm 
the possibility of larger values of D p occurring due to transit-time damping reso- 
nance interaction in the presence of isotropic fast-mode waves in comparison to the 
Alfven waves of the same amplitude (cf. Schlickeiser & Miller 1998). The impor- 
tance of quasi-perpendicular fast-mode waves is stressed for the acceleration of high 
velocity particles. 

Key words: cosmic rays - magnetohydrodynamic turbulence - interstellar medium 
- Fermi acceleration 



1. Introduction 

A number of astronomical objects (extragalactic radio sources, super- 
nova remnants, solar flares) emit radiation with non-thermal spectra. 
These emissions are often connected with the existence of a hot turbu- 
lent magnetized plasma providing conditions for particle acceleration 
by MHD turbulence. It was first shown by Hall & Sturrock (1967) and 
by Kulsrud & Ferrari (1971) that charged particles can be accelerated 
by MHD turbulence having wavelengths long compared to the particle 
gyration radius. For example, the importance of a Fermi-like acceler- 
ation mechanism was considered for extragalactic radio sources (Burn 
1975; De Young 1976; Blandford & Rees 1978; Achterberg 1979; Eilek 
1979). Stochastic acceleration by plasma turbulence is often suggested 
as the acceleration mechanism for energetic particles during impul- 
sive solar flares (see, e.g. Miller & Ramaty 1987; Steinacker & Miller 
1992; Hamilton & Petrosian 1992). It has been proposed that weak 
(5B/B Q <C 1) MHD waves can simultaneously energize ions by cas- 
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cading shear Alfven waves (Miller & Roberts 1995; Miller & Reames 
1996) as well as electrons by the accompanying cascading fast-mode 
waves through a process known as a transit-time acceleration (Miller 
et al. 1996). The MHD waves are attractive candidates because they 
can be produced by either large-scale restructuring of the magnetic 
field, which presumably occurs during the flare energy release phase, 
or by a shear in the plasma bulk velocity (Roberts et al. 1992), which 
is likely to be found in regions of reconnection-driven plasma outflows 
(see, e.g., LaRosa & Moore 1993, Forbes 1996). 

The interaction between the waves and particles is determined pri- 
marily by delta functions, 5(tu — k\\v\\ — nfi), which select those waves 
from a spectrum which are resonant with a given particlefl. At such 
interactions energies vary in a diffusive way, and over a long timescale 
a net energy gain results in the process of stochastic acceleration. For 
n ^ 0, a particle resonates with those waves that are seen at an integral 
multiple of a gyrofrequency and for n = at zero frequency. Then, a 
particle feels a net force which is not averaged out by the phase mixing 
(Stix 1962). The case of n = provides magnetic analogy to Landau 
damping and is associated with a first-order change in 6\B\. This type 
of coupling is not observed in the case of Alfven waves but becomes 
important for the magnetosonic waves containing compressive compo- 
nents of 5B. The interaction between the particle magnetic momentum 
and the parallel gradient of the magnetic field is called transit-time 
damping (cf. discussions by Lee & Volk 1975, Achterberg 1981, Miller 
et al. 1996). 

Recently, Schlickeiser & Miller (1998) presented a quasi-linear 
derivation of cosmic ray transport coefficients in the presence of MHD 
waves, including isotropic fast-mode turbulence. For the isotropic Kol- 
mogorov turbulence they demonstrated that the Fokker-Planck coeffi- 
cients depend both on the transit-time damping and the gyro-resonance 
interactions. For cosmic ray particles with v 3> Va and for vanish- 
ing turbulence cross helicity the momentum diffusion coefficient in the 
fast-mode turbulence is mainly determined by the transit-time damp- 
ing contribution, leading to a more efficient stochastic acceleration in 
comparison to the process in the presence of pure Alfvenic turbulence. 

The aim of the present paper is to study the momentum diffusion 
coefficient D p in the presence of finite amplitude (= 5B ~ 1, here we 
denote 5B/Bq = SB ) magnetosonic and Alfven waves. The influence of 
wave anisotropy and of the wave spectrum slope on D p is considered. 
In the next section we discuss the conditions essential for fast-mode 
wave damping in the plasma. We show that the \ow-(3 plasma can pro- 



The notation is explained in Appendix A 



praca4.tex; 1/02/2008; 5:17; no v.; p. 2 



KLUWER STYLE FILE 



3 



vide conditions with the main damping process being the cosmic ray- 
particles' acceleration. Then, in Section 3, we summarize the results 
of quasi-linear analytic derivations of the momentum diffusion coeffi- 
cient to provide a reference for our numerical modelling. The performed 
Monte Carlo simulations involving derivations of particles trajectories 
in the space filled with finite amplitude fast-mode and Alfven waves are 
described in Section 4. Anisotropic wave distributions are modeled by 
choosing their wave vectors from cones directed along the mean mag- 
netic field B Q . For our simulations we adopt fast-mode or Alfven mode 
turbulence with the flat- (q = 1) or the Kolmogorov (q = 5/3) spec- 
trum within the finite wave vector range (k m i n , k max ). The results are 
presented and discussed in Section 5. We confirm a substantial increase 
of Dp for the fast-mode waves in comparison to the Alfven waves of the 
same amplitude if the nearly perpendicular (k _L B Q ) waves are present. 



2. Dissipation of fast-mode waves in a low (3 plasma 

Fast mode waves can effectively accelerate cosmic ray particles, if other 
dissipation processes are negligible small. We will demonstrate here 
that in gases with low plasma-beta, (3 = 2c 2 s /V^ « 1, this requirement 
is fulfilled. 

The equilibrium intensity of plasma waves results from the com- 
petition of wave generation and wave damping processes. The kinetic 
equation for a 3-dimensional spectral density Wi(k), which denotes the 
wave energy density per unit volume of wavenumber space of the wave 
mode i, is given by the usual conservation equation 

^ -!■%-> . 

where the flux 

m - -D^m ,2,) 

ok 

is expressed as a diffusive term with the diffusion coefficient in 
wavenumber space 

D = k 2 /T s {k) (2.3) 

and the spectral energy transfer time scale T s (k). This approach, 
to describe the evolution of turbulence by a diffusion of energy in 
wavenumber space, was pioneered by Leith (1967) in hydrodynam- 
ics, and subsequently introduced to magnetohydro dynamics by Zhou & 
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Matthaeus (1990). It provides a simple framework to take into account 
- at least approximately - turbulence evolution in space physics appli- 
cations. Zhou & Matthaeus (1990) present a general transport equation 
for the wave spectral density in the case of isotropic turbulence, which 
includes terms for spatial convection and propagation, nonlinear trans- 
fer of energy across the wavenumber spectrum, and a source and sink 
of wave energy. 

For isotropic turbulence one obtains for the associated one- 
dimensional spectral density Wj(fe) = 4irk 2 Wi(k) the simplified dif- 
fusion equation 



dWj d 



dt dk 



4s 



YiWi + Si{k) , (2.4) 




where we included terms for the damping or growth of waves (I^) and 
a wave energy injection and/or sink term Si(k). 

The spectral energy transfer time scale (or the wavenumber diffu- 
sion coefficient (2.3)) depends upon the cascade phenomenology. In the 
Kolmogorov treatment the spectral energy transfer time at a particular 
wavelength A is the eddy turnover time X/5v, where 5v is the velocity 
fluctuation due to the wave. In the so-called Kraichnan treatment the 
transfer time is longer by a factor Va/Sv. Both phenomenologies are 
further discussed in Zhou & Matthaeus (1990) and yield 

,(*•) [ V w > ( Kolm °g° rov ) ; (2 . 5) 

Jjiwi ( Kra i cnnan ) 

where Ub = B 2 /8ir denotes the energy density of the ordered mag- 
netic field. Substituting these transfer time scales into Eq. (2.4), and 
assuming a steady state with no damping, we obtain W{ = Wok~ q , 
where q = 5/3 for the Kolmogorov case and q = 3/2 for the Kraichnan 
phenomenology. The diffusion equation (2.4) in either case is nonlinear. 

Besides the generation of turbulence by kinetic cosmic ray streaming 
instabilities (e.g. Tademaru 1969), wave cascading from low to high 
wavenumbers is very often an important way of producing broadband 
wave spectra. One possibility, discussed e.g. in the context of solar flares 
(Miller & Roberts 1995), is that long- wavelength turbulence results 
from the rearrangement of large-scale magnetic fields and/or a shear 
flow instability, so that it is reasonable to assume the deposition of 
wave energy peaked at long wavelength, probably comparable to the 
physical size of the system, as the primary energy release. Cascading as 
described by the nonlinear diffusion term in Eq. (2.4) will then transfer 



praca4.tex; 1/02/2008; 5:17; no v.; p. 4 



KLUWER STYLE FILE 



5 



this spectral energy to higher wavenumbers, where the waves will be 
able to resonate with progressively lower energy cosmic ray particles, 
until they eventually interact with the charged particles in the tail of 
the background thermal distribution. 

According to Ginzburg (1970) the damping rate of fast magnetosonic 
waves propagating at an angle 9 with respect to the ordered magnetic 
field in a thermal electron-proton background plasma is given by 



T/ | sin 2 6 
16 cos 6 



F(f3,9) 



where 



F(J3,6) 




+ Sexpf-^cos 2 ^ 1 ] 



depends on the plasma beta of the background plasma 



2c| Sim e m p kBT e 



R2 
n 



(2.6) 



(2.7) 



(2.8) 



The function F(/3,9) exhibits an almost perfect "flip-flop" behaviour 

with constant values F{(3 < (3 C ) ~ ^m e /m p ~ 1/43 and F(j3 > (3 C ) ~ 5 

below and above (3 C = 0.5cos~ 2 9 > 0.5. 

The other relevant wave damping process is the transit-time damp- 
ing acceleration of isotropically distributed cosmic rays. According to 
Achterberg (1981) the corresponding damping rate is given by 



cos 9 



1 - 



C 2 COS 2 i 



U B 



(2.9) 



where Ub and U p denote the energy densities of the background mag- 
netic field and the cosmic ray particles, respectively. 

For small wave intensities we can neglect the cascading of waves, 
so that according to Eq. (2.4) the equilibrium fast-mode intensity is 
simply given by 



W f (k) 



_ S f (k) 



(2.10) 



T t + T r 

Calculating following Achterberg (1981) the ratio of the two damping 
rates, 



R = 



V U 



n 

c 2 cos 2 9 



1 -2 



tt^Va 



P 1/2 F(P, 



(2.11) 
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we find that this ratio is independent of wavenumber k, and for small 
plasma beta (3 < 0.5 almost independent from the propagation angle 6 
since Va « c. The ratio is well approximated by 



Values of R < 1 indicate that the fast-mode wave energy is damped by 
accelerating cosmic ray particles. The condition R < 1 translates into 
a condition for the plasma beta 



where T5 = (T/10 5 K) denotes the plasma temperature in units of 10 5 K. 
In the case of equipartition between magnetic field and cosmic rays, 
U p /Ub — 1 we find values of R less unity if (3 < O.Oly/T^. In conclusion, 
in a low /3-plasma the fast-mode waves are predominantly dissipated 
by accelerating cosmic rays by transit-time damping. The inclusion of 
wave cascading does not qualitatively modify this conclusion as the 
numerical solutions of Miller et al. (1996) indicate. 

3. Quasi-linear momentum diffusion coefficient 

The quasi-linear theory treats the effect of the weakly perturbed mag- 
netic field as perturbations of orbits of particles moving in the aver- 
age background field. Schlickeiser (1989) considers quasilinear trans- 
port and acceleration parameters for cosmic ray particles interacting 
resonantly with the Alfven waves propagating along the average mag- 
netic field. The transport equation can be derived from the Fokker- 
Planck equation by a well-known approximate scheme (Jokipii 1966, 
Hasselmann & Wibberenz 1968) which is commonly referred to as the 
diffusion-convection equation for the pitch-angle averaged phase space 
density. The electromagnetic fields generated by MHD waves enter into 
the equation through the Lorentz force term. For fast {v S> Va) cosmic 
ray particles, a vanishing cross helicity state of the Alfven waves and 
the power-law turbulence spectrum with q > 1, 




(2.12) 




(2.13) 



W i {k) = W l {5B i fk- q for (k, 



'mm 




(3.1) 



the momentum diffusion coefficient 




(3.2) 
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where 



Wo = a _ q \ — q for q>l 



and 



Wo= f ln hnaA 1 for q = 1 
V "min J 

Recently, Schlickeiser Sz Miller (1998) considered cosmic ray par- 
ticles interacting with oblique fast-mode waves propagating in a low- 
(3 plasma. In the cold plasma limit the fast and slow magnetosonic 
waves degenerate to the fast-mode waves with the dispersion relation 
u) 2 = V\k 2 . The momentum diffusion coefficient (in the original paper 
'02' for our Dp) and the spatial diffusion coefficient nu can be calculat- 
ed as respective pitch angle averages of the Fokker-Planck coefficient 
-CW(^'P)- Adopting isotropic fast-mode turbulence with a power-law 
turbulence spectrum (3.1) they obtained 



-" 2 > (f)Vr '[frO.) + /G (,)] , (3.3) 

where D w includes a sum of the transit-time damping (fx) and the 
gyroresonance (fa) interaction contributions. A considered form of /t 
admits the pitch angle scattering by transit-time damping of super- 
Alfvenic particles with pitch-angles contained in the range e < \n\ < 1, 
where e = Va/v. In the interval < e, where no transit-time damping 
occurs, the gyroresonance interactions provide a small but finite con- 
tribution to the particle scattering rate. As a result the momentum dif- 
fusion coefficient D p is mainly determined by the transit-time damping 
interactions and the spatial diffusion coefficient by the gyroresonance 
interactions. For 1 < q < 2 fast-mode turbulence spectra (Eq. 3.1) they 
obtained 

2^2 



where 



^^(f)W-^ln^ , (3,) 



Cl ~ 2 (4-^3(1 + 1) • 
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4. Description of simulations 

The approach applied in the present paper for modelling the particle 
momentum diffusion is based on numerical Monte Carlo simulations. 
The general procedure is quite simple: test particles are injected at 
random positions into a turbulent magnetized plasma and their trajec- 
tories are followed by integration of particle equations of motion. Due 
to the presence of waves, particles move diffusively in configuration 
and momentum space . By averaging over a large number of trajecto- 
ries one derives the diffusion coefficients for turbulent wave fields. In 
the simulations we consider relativistic particles with v S> Va and use 
dimensionless units (cf. Appendix A): 6B = SB/B a for magnetic field 
perturbations, l/f2 for time, k/k res for wave vectors and p„Q for the 
momentum diffusion coefficient. Below, in all our simulations we adopt 
Va = lCT 3 c and we consider mono-energetic particles with velocity 
v = 0.99c. This particular choice allows us to compare the effects of 
different types of turbulence on the particle momentum diffusion coef- 
ficient and, due to the large v, to evaluate the role of wave anisotropy 
in the transit-time damping interactions. 

4.1. The Wave Field Models 

In the modelling we consider a superposition of 384 MHD waves prop- 
agating oblique to the average magnetic field B Q = B Q e z . The wave 
propagation angle with respect to B Q is randomly chosen from a uni- 
form distribution within a cone ('wave-cone') along the mean field. For 
a given simulation two symmetric cones are considered centered along 
B Q , with the opening angle 2a, directed parallel and anti-parallel to the 
mean field direction. The same number of waves is selected from each 
cone in order to model the symmetric wave field. Related to the i-th 
wave, the magnetic field fluctuation vector SB^ is given in the form: 

SB^ = SB^ sin(jfeW • r - J l h) , (4.1) 

where 5Bo is a constant vector selected for a given wave. The electric 
field fluctuation related to a particular wave is given as SE^ = — V"W x 
SB^K For Alfven waves ('A') we use the dispersion relation 

u\ = k\Vl , (4.2) 

where Va = B j is the Alfven velocity in the field B Q . The wave 
magnetic field polarization is defined by the formula 

5B A = 5B A (k, lo a ) {k x e z ) fej 1 . (4.3) 
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In a low-/? plasma the fast-mode magnetosonic waves ('M') propagate 
with the Alfven velocity and the respective relations are: 

u 2 M = k 2 V 2 A , (4.4) 

5B M = 6B M (k,uj M ) (k x (k x e z )) fc -1 /cj 1 . (4.5) 

One should be aware of the fact that the considered turbulence model 
is unrealistic at large 5B and the present results can not be considered 
as the exact quantitative ones. In particular, in the presence of a finite 
amplitude turbulence the magnetic field pressure is larger than the 
mean field pressure and the wave phase velocities can be greater than 
Va(B ) assumed here. 

4.2. Spectrum of the turbulence 

In our simulations we consider the power-law turbulence spectrum in 
the wave vector range (k m i n , k max ). The amplitude of the irregular 
component of the magnetic field, obtained from the energy density 
defined in equation (3.1), can be written as 

5B(k) = 6B(k mm )(-^-y q/2 , (4.6) 

where k m i n = 0.08 (k max = 8.0) corresponds to the considered longest 
(shortest) wavelength, respectively, and q is the wave spectral index. In 
the present simulations we consider the flat spectrum with q = 1 and 
the Kolmogorov spectrum with q = 5/3. We included the flat spec- 
trum because of our earlier simulations (Michalek &: Ostrowski 1996), 
where the Alfven waves with q = 1 were considered. It provides also 
a convenient limiting reference for comparison of results for different 
wave spectra flatter than the Kolmogorov one (e.g. the Kraichnan spec- 
trum or the perturbations' spectra at wavelengths longer than the ones 
characteristic for the inertial range of the MHD cascade) . On the oth- 
er hand such a turbulence spectrum is very convenient for numerical 
simulations due to presence of a substantial power in short waves. For 
the flat spectrum the wave vectors are drawn in a random way from 
the respective ranges: 2.0 < k < 8.0 for 'short' waves, 0.4 < k < 2.0 
for 'medium' waves and 0.08 < k < 0.4 for 'long' waves. The respec- 
tive wave amplitudes SB Q are drawn in a random manner so as to keep 
constant 

384 

Y,(SB^?^5B\ (4.7) 
i=i 
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where 5B is a model parameter, and, separately in all mentioned wave- 
vector ranges 

128 256 384 AR 2 

E(^ } ) 2 - E (^) 2 - E (<^) 2 - 4- ■ (4-8) 

i=l i=129 i=257 

Thus the wave energy is uniformly distributed over the considered 
wave-vector ranges. As a second, more realistic turbulence model we 
consider one involving the Kolmogorov spectrum. The observed spectra 
in interplanetary space often have such a form (e.g. Jokipii 1971) and 
this case is most often discussed in the literature. Here all 384 wave 
vectors are drawn in a random manner from the whole considered range 
(0.08 < k < 8.0), but the amplitudes 5B^ are fitted according to the 
Kolmogorov distribution (Eq. 4.6 with q = 5/3) and scaled to keep the 
formula (4.7) valid. In such turbulence most of energy is carried by long 
waves. 

In the discussion below we will consider four different turbulence 
fields labeled as follows: 

i. Alfven waves with the flat spectrum - AF, 

ii. Alfven waves with the Kolmogorov spectrum - AK, 

iii. Fast-mode waves with the flat spectrum - MF, 

iv. Fast-mode waves with the Kolmogorov spectrum - MK. 



5. Results 

In the present numerical simulations we consider a number of physical 
situations - turbulence models with finite wave amplitudes and a chang- 
ing degree of wave anisotropy - not described by analytic means. The 
results are illustrated in figures 1 and 2. First, the simulated momen- 
tum diffusion coefficients, D p , for the Alfven and the magnetosonic 
turbulence are presented in Fig. 1.: in the upper panels variation of 
D p versus the perturbation amplitude 5B for the Alfven turbulence 
and in the middle panels for the fast-mode turbulence. To enable com- 
parison these results are superimposed in bottom panels. The derived 
momentum diffusion coefficients are given for different wave-cone open- 
ing angles. The models with a = 0°, a = 40° and a = 90° represent 
degenerated one-dimensional, anisotropic and isotropic wave vector dis- 
tributions, respectively. For clarity, only the results for the mentioned 
three values are presented, but we also performed simulations for other 
intermediate wave cone openings: a = 30° and 60°. The results of these 
simulations are consistent with the ones described below. 
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At Fig. 1, in all cases a systematic increase of D p with 5B can be 
observed. Usually, it follows the quasi- linear relation D p oc 5B 2 , except 
for the magnetosonic waves with a = 90° , where the increase rate seems 
to be smaller at least at large amplitudes, where it is roughly oc 5B 15 . 

For the Alfven wave turbulence an increase of the wave-cone opening 
angle a providing waves with smaller phase velocities leads to decreas- 
ing Dp. The trend is independent of the considered turbulence ampli- 
tude and the considered spectrum. In the case of fast-mode waves a 
more complicated relation is observed. For the flat turbulence spec- 
trum a small increase of the angle a does not lead to a significant 
variation of D p , and the change is to smaller values of D p , like for the 
Alfven waves. Only the appearance of waves with wave vectors nearly 
perpendicular to the mean magnetic field changes this trend leading 
to an increase of D p (by a factor of three for a = 90° in comparison 
to a = 0° at 5B = 0.2). For the Kolmogorov turbulence the wave 
cone opening-angle increase in the MK model is always followed by an 
increase of D p (in our simulations small values of a, between 0° and 
30°, were not considered) and again it reaches the maximum in the 
presence of perpendicular waves. 

In order to explain this non-monotonic behavior one can refer to 
the quasi-linear derivations of Schlickeiser & Miller (1998). For n = 
the resonance condition for the transit-time damping may be written 
as vu = w/ku and, with the dispersion relation (4.4), Va/vu = k\\/k. 
It is clear that for Va <C v particles can efficiently interact with waves 
at a wide range of vy only if the waves with k± 3> fey are present. 
This fact can explain the effective transit-time damping interactions 
for magnetosonic waves with isotropic distribution. Then, the observed 
Dp increase results mainly from the presence of a small fraction of 
waves propagating quasi-perpendicular to the mean magnetic field. One 
should note that our simulations were performed for relativistic parti- 
cles with v = 0.99c 3> Va, where the effect is quite pronounced. 

In order to verify this behaviour in more detail we considered inter- 
action of such relativistic particles with waves propagating in narrow 
pitch angle ranges with respect to the mean magnetic field: we per- 
formed simulations for wave vector inclinations selected from the angu- 
lar ranges (±5°) around the cones with a = 0°, 45° and 85° (Fig. 2). 
It is clear from the figure that only the presence of waves with wave 
vectors perpendicular to the mean magnetic field leads to a significant 
increase of D p . For these waves the transit-time damping resonance 
provides the wave-particle coupling, while the waves propagating near- 
ly parallel to the average magnetic field contribute preferentially due to 
the gyroresonance. Of course, our reference to the calculations involv- 
ing various resonances is only approximate, as in the presence of high 
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Flat spectum Kolmogorov spectrum 




0.00 0.20 0.40 0.60 0.80 1.00 1.20 0.00 0.20 0.40 0.60 0.80 1.00 1.20 

<5B <5B 



Figure 1. Variation of the momentum diffusion coefficient D p , given in units of 
p 2 o 0. o , versus the perturbation amplitude SB and the wave propagation anisotropy. 
The results are presented for the flat and the Kolmogorov spectra, separately for 
the Alfven wave and the fast-mode turbulence. For comparison the results for the 
Alfven wave turbulence (thin lines) and the fast-mode turbulence (thick lines with 
indicated simulation points) are superimposed at the bottom panels. 
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20000 40000 60000 50000 100000 
Lime 

Figure 2. Examples of the simulated D p , given in units Poilo, versus the simulation 
time for waves selected from the respective (0° ± 5°, 45° ± 5°, 85° ± 5°) pitch angle 
ranges for both types of turbulence spectrum at SB = 0.6. 

amplitude waves the quasi-linear delta form of interaction changes into 
the 'broadened resonance', as discussed by Karimabadi et al. (1992). 
The simulation errors can be evaluated from comparison at Fig-s 1,2 
of the 'M' and 'A' results for a = 0°, which should coincide. 



6. Final remarks 

We considered the momentum diffusion coefficient D p in the presence 
of oblique Alfven and magnetosonic waves with a wide range of ampli- 
tudes, from medium, up to non-linear ones. The influence of the degree 
of wave anisotropy and the waves' spectral index was also studied. As 
expected, in all cases a systematic increase of D p with the wave ampli- 
tude is observed. We confirm a substantial increase of D p for the fast- 
mode waves in comparison to the Alfven waves of the same amplitude, if 
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nearly perpendicular fast-mode waves are present. The effect is caused 
by the transit-time damping interactions which occur in the presence 
of magnetosonic waves containing a compressive component of SB. For 
the isotropic fast-mode Kolmogorov turbulence (cf. Eq-s 3.2,4) D p is 
expected to be on a factor of ln(V/VA), multiplied by a term depend- 
ing on the index q, larger than the one for Alfven waves propagating 
parallel to B Q . The ratio of these values is about three in our simula- 
tions for small 5B, which confirms (within the simulation errors) the 
quasi-linear result. In simulations for the isotropic Kolmogorov spec- 
tra Dp IK /Dp K ~ 10 is roughly constant for a wide range of waves 
amplitudes. For the isotropic MK and MF models one observes that 
the transit-time damping scattering weakens with an increasing 5B. 
For these models, D p oc (SB) 1 ' 5 , when some particles can become time- 
ly trapped between the compressive waves and the acceleration due to 
the transit-time damping becomes less efficient. 

Fast-mode waves can effectively accelerate cosmic ray particles if 
other dissipation processes heating the plasma are negligibly small. In 
Section 2 we demonstrated that for a low beta-plasma, f3 « 1, the fast- 
mode waves can be predominantly dissipated by accelerating cosmic 
rays through the transit-time damping mechanism. In addition, the 
momentum diffusion enhancement due to presence of fast-mode waves 
could also work effectively in a volume with the turbulence generation 
force acting. For example, in the vicinity of a strong shock, where the 
conditions with U p 3> Ub can occur (cf. Drury 1983), or in a region of 
magnetic field reconnection the required fast-mode perpendicular waves 
are expected to be effectively created. One should remember that the 
waves' phase velocities are expected to be larger in realistic turbulence 
than Va(B q ) assumed here. Therefore our results for large 5B should 
rather be considered as the lower limits for D v . 
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Appendix 

A. Summary of notation 

B = Bo + <5B - a magnetic induction vector 

Bo - a regular component of the background magnetic field 

c - the light velocity 

c s - a velocity of sound 

Dp = a2 - a momentum diffusion coefficient 

Dpp- a pitch-angle diffusion coefficient 

e - a particle charge 

E - an electric field vector 

F(k) - an energy flux of waves 

k - a wave vector 

k res = 2n /r g - a resonance wave vector 

k\\ - a wave vector component along B Q 

m e - the electron mass 

m p - the proton mass 

m - a particle mass 

n e - a concentration of electrons 

p - a particle momentum vector (\p\ = p) 

q - a wave spectral index 

r g = c/\Q\ - a particle maximum Larmour radius 
Rl = v/\Q\ - a particle Larmour radius 

R = ryr r 

Si(k) - an injection or a sink term for wave energy 

Ub = Bq/8it - the energy density of the ordered magnetic field 

Up - the energy density of cosmic ray particles 

v - a particle velocity vector (\v\ = v) 

Va ~ the Alfven velocity in the field Bq 

vu - a particle velocity along the mean magnetic field 

Wi(k) - a 3-dimensional spectral density for the wave mode 'i' 

a - a wave-cone opening angle 

(3 = 1(? s jV\ - a plasma-beta parameter 

<5B - a turbulent component of the magnetic field 

Sv - a characteristic velocity fluctuation of the wave 

7 = (1 — v 2 /c 2 ) -1 / 2 - the particle Lorentz factor 

Ti — a damping or a growth term for the wave mode l V 

r r - a damping rate due to the transit-time acceleration 
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r t - a damping rate in a thermal electron-proton plasma 

ijj - a wave frequency 

•jJA - a frequency of the Alfven wave 

ujm - a frequency of the magnetosonic wave 

Q, = eB/^mc - a particle angular velocity 

Q = eB /mc 

li = cos 9 

9 - a wave propagation angle with respect to B Q 
G - a momentum pitch-angle with respect to Bo 
r s (k) - a spectral energy-transfer time scale 
e = V A /v 
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